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Resonances in the superconducting properties, in a regime of crossover from BCS to mixed Bose- 
Fermi superconductivity, are investigated in a two-band superconductor where the chemical potential 
is tuned near the band edge of the second mini-band generated by quantum confinement effects. 
The shape resonances at T — in the superconducting gaps (belonging to the class of Feshbach- 
like resonances) is manifested by interference effects in the superconducting gap at the first large 
Fermi surface when the chemical potential is in the proximity of the band edge of the second 
mini-band. The case of a superlattice of quantum wells is considered and the amplification of the 
superconducting gaps at the 3D-2D Fermi surface topological transition is clearly shown. The results 
are found to be in good agreement with available experimental data on a superlattice of honeycomb 
boron layers intercalated by Al and Mg spacer layers. 
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I. INTRODUCTION 

The main physical properties of the Bardeen-Cooper- 
Schrieffer - Bose-Einstein crossover (BCS-BEC crossover) 
at zero temperature for a system of 3D fermions in- 
teracting via a contact pairing interaction can be de- 
scribed with continuity by the BCS (mean-field) equa- 
tions. Indeed, as shown by Leggett, 1 the ground-state 
BCS wave function corresponds to an ensemble of over- 
lapping Cooper pairs at weak coupling (BCS regime) 
and evolves to molecular (non-overlapping) pairs with 
bosonic character as the pairing strength increases (BEC 
regime). The crucial point is that the BCS equation for 
the superconducting gap has to be coupled to the equa- 
tion that fixes the fermion density: with increasing cou- 
pling (or decreasing density), the chemical potential fi 
results strongly renormalized with respect to the Fermi 
energy Ep of the non interacting system, and approaches 
— \Eb, where Eg is the molecular binding energy of the 
corresponding two-body problem in the vacuum. In re- 
cent years, the experimental realization of the BCS-BEC 
crossover in trapped ultracold fermion gases allowed for a 
quantitative test of the validity of increasingly elaborated 
theoretical approaches: at resonance (the so called "uni- 
tarity limit"), where the lack of a small parameter does 
not permit the systematic control of the approximations, 
the BCS theory at zero temperature gives a gap and 
a chemical potential which are overestimated by about 
30% as compared with the results of QMC calculations 
and the experimental findings, while already the inclusion 
of pair fluctuations at the (non-self-consistent) t-matrix 
level improves the comparison for the gap and /i, with de- 
viations not larger than few percents, as shown by quan- 
titative comparison between theoretical predictions and 
experimental results for the BCS-BEC crossover. 2 The 



mean field approach at zero temperature, within the local 
density approximation, 3 or in the form of Bogoliubov-de 
Gennes equations, 4 has been widely used in the context of 
ultracold atoms to study nonuniform configurations both 
in homogeneous space and in the confining potential of 
the trap, as, e.g., quantized vortices. On the basis of 
the wide range of quite successful applications discussed 
above, in this paper we shall apply the BCS approach 
(with the inclusion of the density equation) to the case 
of multiband superconductivity with the chemical poten- 
tial tuned near a band edge. This approach, which has 
the advantage of being numerically feasible, is also apt to 
give a good overall insight into the resonant and crossover 
phenomena which are the object of the present work. 

Multigap superconductivity in a system with many 
mini-bands of different symmetry has been studied in 
electronic systems with relevant quantum size effects: 
single slab, 5,6 single tube, 7 single wire, 8,9 or single quan- 
tum dots. 10 The electronic structure generated by the 
quantum confinement effects gives rise to electron wave- 
functions at the Fermi level near a band edge that are 
strongly affected by the detail of the quantum confine- 
ment. Recently, experimental evidence has been accu- 
mulated for superconductivity modulated by quantum 
confinement effects in a single quantum well and in a 
single quantum dot. 11-15 

Also superlattices of quantum wires or quantum wells 
show multiband superconductivity near a band edge. 16-20 
The three dimensional (3D) structure of the superlat- 
tice has the advantage to suppress quantum fluctuations, 
that reduce the critical temperature in low dimensions, 
while keeping the key features of quantum confinement 
effects for superconductivity. While there are works on 
the theory of shape resonances in two dimensional (2D) 
superlattices of quantum wires, 21 ~ 23 called "superstripes" 



in cuprate superconductors, 24-26 the theoretical investi- 
gation of shape resonances in superlattices of quantum 
wells is missing. These are now of high interest since a 
3D superlattice of quantum wells provide the simplest 
case of a 3D system showing multiband superconduc- 
tivity near a band edge. Moreover, it has been rec- 
ognized that diborides 27 , intercalated graphite 28-30 and 
pnictides 31 ' 32 are practical realizations of a superlattice of 
superconducting layers at atomic limit, where interband 
pairing is an essential ingredient for high-temperature 
superconductivity 16-20 

In this paper, we study the BCS-BEC crossover which 
occurs in a two-band system, when the chemical poten- 
tial jjl is tuned in a narrow energy range around the edge 
E e dge of the second mini-band. This description applies 
to the situation when the first two mini-bands produced 
by quantum size effects in a superlattice of quantum wells 
are well separated, i.e., when the electron hopping be- 
tween wells is small enough, to that the corresponding 
transversal band dispersion £ is smaller than the energy 
separation between the mini-bands. 

In Fig. 1 we show that a new 3D Fermi surface (FS) 
opens when /x is increased above E e d ge , and the electron 
gas in the metallic phase undergoes an electronic topo- 
logical transition (ETT), also called a Lifshitz quantum 
phase transition, of the type "appearing or disappearing 
of a new Fermi surface spot". 19 ' 20,33 When fj, reaches a 
higher energy threshold, -E3D-2D, the electronic struc- 
ture of the superlattice undergoes a second ETT, the 
3D-2D ETT where one FS changes topology from 3D to 
2D (i.e., from "spherical" to "cylindrical") or vice versa, 
called "the opening or closing of a neck in a tubular 

FS v 19,20,27,33 Tmg ETT ^ typical Q f j) stacks Q f meta l_ 

lie layers, ii) multilayers, or iii) superlattices of quantum 
wells, and therefore it is a common feature of all exist- 
ing high-temperature superconductors and novel materi- 
als synthesized by material design in the search for room 
temperature superconductivity. In practical cases, the 
tuning of /1 at an ETT can be controlled by means of: 
i) pressure; ii) the misfit strain 34-36 between the super- 
conducting layers and the spacers; iii) the ordering of 
dopants in the spacers 26 with related striped superstruc- 
tures; iv) the effective charge density, via gate voltage 



methods; 



.37-39 



the charge transfer between the super- 



conducting layers and the spacers; vi) the thickness of 
the spacers. 

As stated above, we focus on the superconducting 
properties of a system in which the chemical potential 
is varied in a narrow energy range around the band edge, 
within a simple model for a two-band superconductor 
with a first and a second band of different (even and 
odd) parity and different topology, which can be varied 
by tuning the chemical potential. This system provides 
a particular case of modulated topological multiband su- 
perconductor with multiple different winding numbers. 
We consider the set of coupled BCS equations to study 
the evolution of the properties of this two-band system 
of fermions with contact pairing interactions and a pair- 




FIG. 1. Pictorial view of the evolution of the Fermi sur- 
faces (FS) for the superconductor made of non-hybridized 
two-bands by moving the chemical potential from panel (a) to 
panel (d) so that it crosses two Lifshitz electronic topological 
transitions (ETT). The first ETT occurs moving Ef across 
the band edge E e d ge of the second band so that the supercon- 
ductivity goes from the single FS in panel (a) with a single 
condensate to the two FS in panel (b) with two condensates: 
the first (1) has a 2D topology and the second (2) has a 3D 
topology. Here, £ is the transversal band dispersion. Chang- 
ing Ef the system crosses the critical energy E3D-2D where 
the second FS undergoes a 3D-2D ETT shown in panel (b) 
changing its topology: the second closed 3D FS [panel (b)] 
becomes the tubular 2D FS in panel (d). The first large 2D 
FS (1) remains nearly constant when the chemical potential 
is moved. 



ing energy cutoff ljq on the order of magnitude of the 
distance between jjl and the bottom of the second band. 
Moreover, resonant and crossover phenomena which are 
the object of this work, are amplified when the transver- 
sal dispersion £ is of the same order of magnitude of cjo- 
Therefore, in the following discussion, we concentrate our 
analysis on systems that are engineered by material sci- 
ence techniques in such a way that £ = ojq. 

The condensate in the first large 2D FS, in the stan- 
dard BCS approximation (Ep/ujo 3> 1, i.e., Ep far 
from band edges), coexists with a second small FS, 
where the standard approximation breaks down [(Ep — 
E e dge )/ UJ o "C 1, i.e., Ep close to the band edge where a 
new FS disappears or appears, or near the 3D-2D ETT, 
at the closing or opening of a neck in a tubular FS 33 ]. 

When the chemical potential is located near the band 
edge, it may be driven below the bottom of the con- 
duction band by superconducting correlations, in a way 



similar to the case of semimetals. 41 The phenomenon of 
quantum shape resonances 16-20 occurs when the chemi- 
cal potential falls within an energy range wo near an ETT 
point. The shape resonances are due to the configuration 
interaction between paired states in the large FS, with a 
large group velocity, and quasi stationary states near the 
ETT in the second small FS. The shape resonances in 
superconducting gaps near an ETT are analogous to the 
scattering resonances due to configuration interaction be- 
tween a closed and an open channel described in nuclear 
physics by Majorana and Feshbach, and in atomic physics 
by Fano, as reported in a recent work. 42 

In our system, we capture the crucial ingredients of 
the interplay of the two bands and of the two set of 
wave functions (which will be shown to be responsible for 
the shape resonances in the gap and in the critical tem- 
perature), without facing the complications beyond BCS 
theories (as the Eliashberg theory, or corrections beyond 
Migdal theorem, etc.) which are considered in more con- 
ventional (phonon- or magnon-mediated) strong-coupling 
superconductors, or the t-matrix corrections. These lat- 
ter would be the dominant corrections in the regime of 
band-edge crossing, in particular for the evaluation of the 
critical temperature, but should not qualitatively change 
our results. Indeed, the Migdal-Eliashberg approach to 
evaluate the superconducting properties taking into ac- 
count the frequency dependence of the two-particle in- 
teraction (retardation effects) can be applied only when 
the Fermi energy is much larger than the typical en- 
ergy of the interaction, i.e., when the Migdal theorem 
holds. In multiband superconductors, when approaching 
a band bottom, the two energy scales become compara- 
ble and the Migdal theorem cannot be applied (a small 
parameter controlling the diagrammatic expansion does 
not exist anymore) , requiring the inclusion of much more 
complicated diagrams, as vertex corrections, in determin- 
ing the self-consistent equations for the superconducting 
properties, which is beyond the scope of our work. 

We point out that, in this paper, we shall not deal with 
the problem of superconducting fluctuations in multi- 
band systems, and refer to the different regimes depend- 
ing on the position of the Fermi energy only. Moreover, 
we recall that a pure BEC limit with the condensation of 
point-like bosons, cannot be achieved when a pairing in- 
teraction with a finite energy cutoff is considered in the 
effective Hamiltonian. When the Fermi energy crosses 
both the band edge, the condensate in the newly disap- 
pearing or appearing FS undergoes a crossover from a 
mixed Bose-Fermi regime to the BCS regime. Indeed, 
in the range — 1 < (Ep — E e d ge )/^o < 0, the electron 
states associated with the newly appearing FS are un- 
occupied in the normal state, and, in the presence of a 
strong enough pairing interaction, a condensate of boson- 
like tightly bound pairs is formed below T c . However, we 
call this regime the mixed Bose-Fermi regime, because 
the tightly bound pairs pairs are never genuinely boson- 
like, due to the presence of the large band. A preliminary 
study of the superconducting fluctuations 43 reveal that 



these have a mixed character, propagating (like in the 
BEC regime) and damped (like in the BCS regime). 

The first unconventional BCS regime occurs in the 
range < (Ep — E e d ge )/^o < 1, where all the few elec- 
trons in the newly appearing mini-band, condense below 
T c . Indeed, in the case under discussion, £ = loq, the stan- 
dard approximation (Ep — E e d ge )/ujQ ^$> 1 breaks down. 
In this range, the new appearing FS has a 3D topology 
while the large FS has a 2D topology. 

A second unconventional regime occurs around the 3D- 
2D ETT where the superconductivity phase in the range 
< (Ep — E e dge )/wo < 2 arises because of the config- 
uration interaction between pairing channels in different 
bands with different condensate symmetry 2D and 3D 
(i.e., with different winding numbers) above and below 
the 3D-2D ETT. 

Finally the system is in the standard multiband BCS 
regime for (E F - E edge )/u > 2. 45 ~ 47 

There is an analogy between the crossover case stud- 
ied here and the BEC-BCS crossover in ultracold Fermi 
gases. While in the ultracold gases the tuning of the en- 
ergy of the bound state of the diatomic molecule above 
or below the continuum is performed by using an ex- 
ternal magnetic field, here, the ETT in the second nar- 
row mini-band can be tuned for example by gate voltage 
techniques, pressure or misfit-strain, inducing shape res- 
onances in the superconducting properties. 

Here, starting from the standard two-band BCS su- 
perconductivity in the clean limit 45 " 47 far from ETT's, 
we focus on the case when the chemical potential falls 
within a narrow energy range near a band edge where 
the Feshbach-like shape resonance in superconducting 
gaps takes place. 19 The shape quantum resonance in 
the exchange-like interband interactionis is generically 
driven by the repulsive Coulomb interaction and in a two- 
band model leads to condensate wavefunctions with op- 
posite signs (a classical quantum phenomenon, recently 
becoming popular in pnictides with the name of s± 
pairing). 19 ' 20,49 ' 50 We show that the control of the ratio 
between the intensity of exchange-like interband pairing 
and intraband Cooper pairing, by material design tech- 
niques, is crucial. We discuss the case of different attrac- 
tive intraband coupling strengths in the first and second 
mini-band and we determine the evolution of the super- 
conducting gaps by changing the interband exchange-like 
pairing in the crossover regime from the BCS regime, well 
above the band edge to the mixed Bose-Fermi regime, be- 
low the band edge. 

We find that the direct evidence for the quantum in- 
terference effects between pairing channels is provided by 
minima in the gap parameter for electrons in the large FS 
and we show the shift to these minima by changing the 
interband pairing strength. We report evidence for two 
different regimes where in the first the interband pairing 
dominates while in the second the intraband dominates. 
These regimes are separated by a critical value of the 
interband pairing where the gap ratio becomes indepen- 
dent on the variation of the chemical potential. Finally, 



we are able to compare the calculations with the avail- 
able experimental data for the evolution of the gaps near 
the band edge of a 2D FS in doped diborides, where in- 
traband dominates, but interband pairing is essential for 
understanding the evolution of the gaps with the (doping 
dependent) critical temperature. The comparison with 
available experimental data for doped diborides in a wide 
range of doping provides a good test of the present the- 
ory. 



II. THE SUPERCONDUCTING GAPS 

In a multiband system, near the band edge in the £- 
th band, where Vk-E^ = 0, the energy of an electron 
can be approximated by a free-electron dispersion law 
Ei.k — Ei + (k 2 /2mi), where mi is the electron effective 
mass at the band edge. However, within our model for a 
supcrlattice of quantum wells disposed along the z axis, 
while this approximation is valid to describe the electron 
dispersion in the x, y plane of the superconducting lay- 
ers, it is certainly not valid in the z direction, when the 
transversal band dispersion £ is on the order of energy 
cut off of the pairing interaction luq. In this case, an 
anisotropic band with weak dispersion in the z direction 
and larger dispersion in the x, y plane above the band 
edge of the second band should be adopted, 
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where E{k z ) is the actual energy dispersion in the peri- 
odic potential of the superlattice and mn is the effective 
mass in the x, y plane. This situation is obtained within a 
model where a free electron gas is confined in a potential 
which is periodic in the z direction, 
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W(z)= J2 Wb(z-nd), 
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where W b (z) = -V b for \z\ < L/2 and W b (z) = for 
L/2 < \z\ < d/2, L is the width of the confining well 
and d is the periodicity of the superlattice in the z direc- 
tion. This periodic potential mimics the phenomenology, 
e.g., of the pnictides, diborides, and stacks of graphenc 
layers made of a superlattice of stacked planes. The con- 
fining potential generates a band structure organized in 
mini-bands. Each mini-band has a 3D character, with 
closed isoenergetic surfaces, near the lower band edge 
E e dge = E^l, and a 2D character, with isoenergetic sur- 
faces open in the z direction, above some energy thresh- 
old Ezd-2D = Ei.t- The present model is apt to describe 
quantum interference phenomena between different scat- 
tering channels in a large and a small FS which are the 
object of this work, when \x is tuned near the bottom of a 
f-th mini-band, within a window of width 4wo- In the en- 
ergy window of width 2wo the first band, with a large 2D 
tubular FS, has a constant density of states (DOS) N\. 



The second FS disappears (or appears) as Ep crosses 
the level E e d ge = Ei^. The second FS changes from 
tubular 2D to closed 3D topology as Ep crosses the level 
E3D-2D = Ei y T, as it is shown in Fig. 1. Near the edge, 
the DOS of the second FS, N2, has the typical 3D be- 
havior shown in Fig. 2a in the corresponding range of 
charge density in the metallic layers between p e d ge and 

P3D-2D = Pedge + Pc- 

The crossover from 2D to 3D can be described in our 
model by changing from an infinite potential barrier be- 
tween the planes of the superlattice in the z direction, 
provided by spacer layers, to a finite barrier (Vb). This 
yields a finite hopping term that broadens the sharp dis- 
continuity of the DOS of a pure 2D band. This broaden- 
ing increases the width of the shape resonance. Moreover, 
it is possible to design artificial superlattice heterostruc- 
tures at optimum shape resonance condition, i.e., where 
the value of the potential barrier V b and its width are 
such that the mini-band dispersion in z direction £ is of 
the order of energy cutoff wo of the interaction, as we 
shall assume in the following. 

It must be pointed out that, differently from the case 
of a single 2D slab, 5 in our superlattice, 21 at finite V b a 
genuinely 3D condensate is formed, reducing the effect 
of fluctuations of the superconducting order parameter, 
which suppress the mean-field value of T c . Moreover, co- 
existence of small and large Cooper pairs due to different 
couplings in the two bands is expected to induce an effec- 
tive screening of the strong superconducting fluctuations 
associated with the small pairs. 48 

The pairing interaction is assumed to be originated 
from an electron-electron contact interaction. Once the 
matrix elements between exact eigenstates of the super- 
lattice, V k ' k , , are calculated they turn out to depend on 
the wave vector k z in the superlattice direction. This in- 
duces a structure in the k-dependent interband coupling 
interaction for the electrons that determines the quan- 
tum interference between electron pairs wave functions 
in different mini-bands of the superlattice. 19,20 The pair- 
ing interaction is then cut off at a characteristic energy 

Vui = "0(^0 - |&k|)0(u;o - |&',k'|) (3) 

where k = k z (k' — \a' z ) is the superlattice wavevector, in 
the z direction, perpendicular to the planes, of the initial 
(final) state in the pairing process, and 
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where Nq(Ep) is the DOS at Ep for a free electron 3D 
system, V%d is the volume of the system, 
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and ipijs.(z) are the eigenfunctions in the superlattice of 
quantum wells, normalized so that f d dz\ip£^(z)\ 2 = 1. 
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FIG. 2. (Color online). The crossover regime reached by 
decreasing the charge density in the superlattice from the two- 
band BCS superconductor to the mixed Bose- Fermi regime at 
the lover band edge E2,l, as a function of the electron charge 
density (p - p e dge)/Pc where p edge (p edge + p c ) is the charge 
density where Ef is tuned at E e dge {E-$d-2d)- Panel (a). 
The ratio of the DOS of the second mini-band N2 and the 
DOS of the first mini-band Ni. Upper part of panel (b): The 
relative variation of the Lifshitz parameter (p n — E2,l)/£ as a 
function of electron charge; lower part panel (b): the variation 
of the chemical potential, (p s —p n ) / p n , as a function of charge 
density p, where p n and p 3 are the chemical potentials in the 
normal and superconducting phase, E^,l is the lower band 
edge energy of the second mini-band. 
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FIG. 3. (Color online). Density histograms of the pairing 
interaction matrix elements 7 k ' k / between the first and the 
second mini-bands, for the case with mini-band dispersion 
£ = (jJo- Histograms 11 (blue), 22 (light blue), 21 (violet) 
and 12 (red) shows respectively the pairing interaction matrix 



The use of single cutoff in two-band superconductors has 
been justified in detail by Entel. 51 Here, we assume that 
pairing in our confined geometry is provided by a contact 
interaction with a given characteristic energy range. The 
dimensionless coupling constants Ce t i> = (2Si t £> — l)c°£, 
are positive for £ = £' (intraband Cooper pairing) and 
negative for I 7^ £' (repulsive exchange-like interband 
pairing, with ctjy = ci'.e) and measure the relative in- 
tensity of intraband and interband pairing. 

Here a comment on the effective interaction adopted in 
the present work is in order. Our pairing interaction is 
semi-phenomenological in the sense that we do not derive 
the microscopic expression of the contact pairing interac- 
tion (its strength and energy range being used as fitting 
parameters), but the matrix elements are calculated be- 
tween eigenstates of the actual confining potential, and 
in this sense they are microscopic. The choice of luq on 
the order of the mini-band width £ is justified because 
this corresponds to the optimum amplification case, as 
above discussed. 

In order to determine the gaps and p self-consistently 
at zero temperature we use iterative solving methods for 
the BCS-like equations 
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and 7, 



where M is the total number of wave- vectors k' and p is 
the electron density, starting with initial constant gaps 
and an initial p equal to the value of Ep in the normal 
state. We considering convergence to have occurred when 
the relative variation of the gap and charge density p 
less than 10 ~ 6 . We point out that the gap functions 
A^k depend on the superlattice wave- vector k as well as 
on the mini-band index £. The BCS-like equations self- 
consistently couple the gap function at a given point of 
k-space with the values of the gap function in the entire 
k-space. In the following we discuss the properties of the 
average values of the gap in k-space, unless otherwise 
specified. The above equations can be easily generalized 
to finite temperatures T, but in the remaining part of the 
present section we focus on the case T = 0. 

In the standard BCS theory, where Ep is far from the 
band edges, the relative variation of p going from the 
normal to the superconducting state is expected to be 
negligible. This is not true when p is tuned near the 
band edge of the second band. In fact, our calculation 
yields a significant variation of p in the superconducting 
phase, as a function of the charge density (Fig. 2b). 
A relative variation of p going from the normal to the 
superconducting phase, as large as 10~ 3 , is obtained near 
the band edge and at the 3D-2D ETT, within a range 
4wq. The variation starts to be large, as compared with 
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FIG. 4. (Color online). Calculation of the gap anisotropy, 
i.e., the dependence of the gap in the first mini-band, with 
C22/C11 = 2.17, en = 0.22 and cia/cn = 0.43 as function of 
the wavevector k z in the transversal direction. We show in 
the figure two characteristics cases with different values of the 
averaged gap obtained solving the BCS-like equations at finite 
temperatures T = 5.7K (solid line) and T = 27. 2K (dashed 
line) in the first mini-band. 



the standard BCS result, in proximity of the crossover 
regime below the band edge up to well beyond the 3D- 
2D ETT. 

The Feshbach-like shape resonance regime occurs in 
correspondence of the large variation of \x between the 
normal and superconducting phases shown in Fig. 2b. 
In Fig. 3 we report the distribution of values of the ma- 
trix elements / kk / in the first two mini-bands, for the 
case with mini-band dispersion £ = uj$. The intraband 
distributions of the two bands show different shapes and 
widths and have different range of values. The resulting 
matrix of coupling constants is obtained exclusively from 
the eigenfunctions of the superlattice and is noticeably 
asymmetric. 

The gap parameters in the two bands have a sizable 
dependence on the wave- vector k z as shown in Fig. 4 for 
the gap in the first mini-band, taking the average value of 
the gap over the wave-vectors as a reference value. The 
wave-vector dependence of the gaps is induced by the 
effective interaction V k ' k , entering the BCS equations, 
which by itself has a non-trivial dependence on the wave- 
vectors of the two scattered electrons forming the Cooper 
pairs, as discussed above. The largest values of the gaps 
are found in our calculations for large transversal wave- 
vectors, suggesting that the effective pairing interaction 
is mostly affected by the superlattice at small distances 
of the order of the modulation of the confining periodic 
potential. Therefore, all the superconducting properties 
(coherence length, critical temperature, gaps to T c ra- 
tios) are expected to be influenced by superlattice effects 
via the wave- vector dependence of the gap parameters in 
different bands. 

The variation of the gaps at T = as function of the 
chemical potential in different coupling regimes are plot- 
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FIG. 5. (Color online). The evolution of the gap, averaged in 
the k z direction, in the first mini-band [panel (a)] and in the 
second one [panel (b)] as a function of the Lifshitz parameter, 
with C22/C11 = 0.45 and en = 0.22. The ratio of the gaps as 
a function of the Lifshitz parameter is reported in panel (c). 
Several interband coupling ratio cases are reported: curve 
1 (solid line) shows the case for C12/C11 = —2.73, curve 2 
(dashed line) shows the case for C12/C11 = —1.59, curve 3 
(dotted line) shows the case for C12/C11 = —1.04, curve 4 (dot- 
dashed line) shows the case for C12/C11 = —0.91 and curve 5 
(solid line with open squares) shows the case for C12/C11 = 
—0.68. The curve 3 (dotted line) shows the critical case of 
interband pairing _ c g ltlc<l! — y/cnew/ 0.656 where the gap 
ratio remains constant going from the two-gap BCS Fermi 
case to the mixed Bose-Fermi regime at the band edge. 



ted in Figs. 5 and 6. In both figures, the average value 
of the two gaps taken over the transversal wave-vector 
k z is reported. In Fig. 5 we consider the case where the 
intraband coupling in the second mini-band is smaller 
than the intraband coupling in the first mini-band, i.e., 
C22/C11 = 0.45. In particular, we report the gaps in the 
first [panel (a)] and second [panel (b)] mini-band as a 
function of the chemical potential tuned in the vicinity 
of the bottom of the second mini-band. The ratio be- 
tween the two gaps is reported in the same figure, panel 
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FIG. 6. (Color online). The evolution of the gap, averaged in 
the k z direction, in the first mini-band [panel (a)] and in the 
second one [panel (b)] as a function of the Lifshitz parameter, 
with C22/C11 = 2.17 and en = 0.22. The ratio of the gaps as 
a function of the Lifshitz parameter is reported in panel (c). 
Several interband coupling ratio cases are reported: curve 
1 (solid line) shows the case for cia/cn = —6.09, curve 2 
(dashed line) shows the case for C12/C11 = —4.34, curve 3 
(dotted line) shows the case for C12/C11 = —2.24, curve 4 (dot- 
dashed line) shows the case for C12/C11 = —1.08 and curve 5 
(solid line with open squares) shows the case for C12/C11 = 
—0.43. The curve 3 (dotted line) shows the critical case of 
interband pairing - c c {J Ucal = ^ciic 22 /0.656 where the gap 
ratio remains constant from the two-gap BCS Fermi case to 
the mixed Bose-Fermi regime at the band edge. 



(c). Several interband couplings are considered, span- 
ning from (ci2 = C21 = —0.15) to {c\i = C21 = —0.60). 
Interestingly, a pronounced minimum is present for the 
average value of the gap at the FS of the first mini-band, 
Ai, when the chemical potential approaches the bottom 
of the second mini-band from below. The partial DOS 
and the intraband coupling in the first mini-band do not 
change moving the chemical potential therefore the pres- 
ence of a pronounced minimum supports the existence of 
antiresonances in the superconducting gaps induced by 
the interplay between the wave- vector and band index 
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FIG. 7. (Color online). Panel (a). The variation of the 
critical chemical potential (solid line with filled circles) as 
a function of intraband coupling ratio caa/cii. We observe 
that for the intraband coupling ratio approaching unity the 
critical chemical potential goes to infinity. Panel (b). The 
variation of the critical interband pairing as a function of 
y/ciic 2 2- We observe a linear fit y/cHoE = -0.656 cf"*""' 



(solid line) of calculated Ci" lca points for the cases of y/cuC22 
(open circles) that divide the plane space in two regions: 
the upper one, for y / CiiC22~ > _ c g li * ca! j i n weak inter- 
band regime, where 9(— A2/Ai)/dfi > and the lower one, 
for ^/cnC22 < —Ci2' t ' cal , in strong interband regime, where 
9(-A 2 /Ai)/<9,u < 0. 



dependence of the effective interaction (in which the su- 
perposition of single-particle wave-functions gives rise to 
interference effects) . The depth of the minimum and its 
energy position depend on the interband exchange-like 
pairing. It generally appears below the band edge, where 
the DOS of the second mini-band changes abruptly. It 
occurs well below the band edge for strong interband 
pairing and moves toward the band edge by decreasing 
the interband pairing. The details of the antiresonancc 
depend on the interband coupling, but its presence is a 
robust feature. 

The average value of the superconducting gap at the 
FS of the second mini-band, A 2 , appears below the en- 
ergy of the band edge where the chemical potential enters 
in the second mini-band. It displays a strong interband 
coupling dependence [see Fig. 5, panel (b)]. The ratio 
— A2/A1 reaches a constant energy independent value at 



a critical value c1" tlcal of the repulsive exchange-like in- 
tcrband pairing. The average value of the superconduct- 
ing gap at the FS of the second mini-band, A 2 , shows 



an antiresonancc minimum for | C12 1 < 



^critical 



The 



gap ratio — A2/A1 shows a divergent enhancement mov- 
ing the chemical potential from the two-gap BCS regime 
well above the band edge toward the mixed Bose-Fermi 
regime near the edge for | C12 1 > \c c { 2 %t%cal \- The gap ratio 
— A2/A1 decreases moving the chemical potential from 
the two-gap BCS regime well above the band edge to- 
ward the mixed Bose-Fermi regime reaching a minimum 
at the band edge and increasing below the band edge for 
I c 12 1 < \cl2 lUcal \, as shown in Fig. 5, panel (c). When 
the chemical potential is well inside the second mini- 
band, the two gaps recover the standard behavior of a 
two-band superconductor, with the gap ratio related to 
the coupling ratio C22/C11. Therefore — A2/A1 remains 
smaller than one in the present case. The results show 
clearly that in the crossover regime there is a critical in- 
terband pairing c c™* lcai that separates a first regime of 
strong interband, where the gap ratio — A2/A1 increases 
with moving the chemical potential toward the band edge 
from a second regime of weak interband, where the gap 
ratio — A2/A1 decreases approaching the band edge. 

In Fig. 6 we consider the case where the intraband 
coupling in the second mini-band is larger than the intra- 
band coupling in the first mini-band, i.e., C22/C11 = 2.17. 
In particular, the gaps in the first [panel (a)] and sec- 
ond [panel (b)] mini-band as a function of the chemi- 
cal potential tuned at the bottom of the second mini- 
band are reported, together with the ratio between the 
two gaps in panel (c). The calculations have been car- 
ried out for several interband couplings, spanning from 
weak (ci2 = C21 = —0.51) to strong interband exchange- 
like repulsive coupling (= —1.4). The pronounced an- 
tiresonancc minimum is clearly present for Ai. The ra- 
tio — A2/A1 reported in panel (c) of Fig. 6 reaches a 
constant energy independent value at the critical value 
c cnticai Q £ ^ ne j-epulsive exchange-like interband pairing. 

Our results for the case Onjcw > 1 show that the ra- 
tio — A2/A1 display another relevant feature: at a critical 
position of the chemical potential the ratio between the 
two gaps becomes universal, i.e., independent of the value 
of the interband coupling. The critical chemical potential 
in Fig. 6 is above the bottom of the second mini-band, at 
a distance about ^wo from it. This is again a signature 
of interference effects in the effective pairing, resulting in 
a fixed point between an antiresonance and a resonance 
in the superconducting gaps. In Fig. 6 we observe that 
well above the band edge, in the standard two-gap BCS 
regime, the gap ratio — A2/A1 is larger than one, when 
such is the ratio between the intraband couplings. In- 
creasing the interband pairing we observe that there is a 
value of the chemical potential at which there is a reversal 
of the ratio — A2/A1. 

The critical point of the chemical potential where the 
gap ratio is independent of the interband pairing for 
C22/C11 > 1 depends on the ratio between the intraband 



attractive couplings in the two bands. The variation of 
the critical chemical potential, present only for the case 
C22 > en, depends on the ratio between the intraband 
couplings C22/C11 > 1, as it is shown in Fig. 7 [panel (a), 
dashed blue line] . The critical chemical potential goes to 
infinity when the couplings are equal and decreases to- 
ward the band edge with increasing the ratio between the 
coupling in the second band on the coupling in the first 
band. The critical interband pairing c^ t% increases 
linearly with increasing ^/cnC22, as it is shown in Fig. 7 
[panel (b), solid red line]. The data can be fitted with the 
linear relation y^cncii — — 0.656 Ci" tl , as it is shown 
in Fig. 7 [panel (b)]. 



III. THE CRITICAL TEMPERATURE 

The superconducting critical temperature T c is calcu- 
lated by iteratively solving the linearized gap equation 



A# k = - 
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U'k' 



(8) 



until the non-trivial solution (having eigenvalue 1) is 
reached with increasing temperature. 

Below, we present numerical results for the solution 
of the self-consistent linearized BCS equations which de- 
termines the value of the critical temperature and of the 
chemical potential. We analyze the behavior of the ratios 
2Ai/T c and 2A2/T C , where, according to our notation, 
Ai and A2 are the average values of the wave- vector de- 
pendent gaps on the corresponding branches of the FS at 
T = 0. Note that the gaps to T c ratios are physical quan- 
tities which are very suitable for quantitative comparison 
with experiments. Indeed, gaps at T — (or temperature 
below 0.4 T c in real experimental conditions), and T c can 
be extracted in a single experimental run when measur- 
ing the temperature dependence of the superconducting 
gaps in multiband superconductors. 

We focus here on the case characterized by a coupling 
in the second FS larger than in the first FS. As an exam- 
ple of this regime, we calculate the DOS for the second 
mini-band, T c , 2Ai/T c and — 2A2/T C as functions of the 
reduced Lifshitz parameter (/i— .E2,l)/<*>o, with intraband 
coupling terms fixed to en = 0.23, c-iijcw = 2.17, for 
several values of the interband coupling ratio, as shown 
in the panels (a-d) of Fig. 8. 

The minima of T c , due to the shape antiresonances, 
are all below the band edge energy of the second mini- 
band [panel (a)] while the maximum value of T c , due to 
the shape resonance, is reached near the 3D-2D ETT of 
the second mini-band, on its 2D side [panel (b)]. We ob- 
serve that the shape resonances in the gap to T c ratios 
show clear evidence for the typical asymmetric line shape 
of Fano-Feshbach quantum resonances driven by configu- 
ration interaction between different scattering channels. 
The antiresonance due to negative interference effects ap- 
pears in the range between (p — E2,l)/ uj q = — 1 and 
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FIG. 8. (Color online). The case of strong coupling in the 
second band (the so-called diboride case) with C22/C11 = 2.17 
and en = 0.22. The variable ratio u — N2/N1 (solid black 
line) is shown in panel (a). The critical temperature T c 
[panels (a) and (b)], the ratio 2Ai/T c [panel (c)], the ratio 
— 2/S.2/T c [panel (d)], as functions of the reduced Lifshitz pa- 
rameter (/1 — E2,l)/wo. Several values of the interband cou- 
pling ratio are reported: curve 1 (solid line) shows the case 
for C12/C11 = —6.09, curve 2 (dashed line) shows the case 
for C12/C11 = —4.34, curve 3 (dotted line) shows the case 
forci2/cu = —2.24, curve 4 (dot-dashed line) shows the case 
for C12/C11 = —1.08 and curve 5 (solid line with open squares) 
shows the case for C12/C11 = —0.43. For the largest value of 
C12 considered here, the critical temperature reaches 300K, 
in the range 1< (/x — E2 1 l)/wo <2. At the antiresonance for 
(p — E2,l)/oJo = 0, T c and 2Ai/T c reach the minimum value 
for weak interband interaction. On the contrary, for large in- 
terband interaction, 2Ai/T c and T c reach the minimum value 
at (fi — E2,l)/oJq = — 1, where T c reaches room temperature 
for the chosen set of parameters. 



(/i — -E^.lV^o = 0, where T c has a minimum in panel 
(a) followed by the T c resonance maximum in the range 
(/x- E 2: l)/oj = 1.5. 

For strong exchange-like interband pairing, C12 > 
we find that the ratio 2Ai/T c shows two minima 
and the ratio 2A2/T C shows an increasing divergence in 
the mixed Bose-Fermi regime. Decreasing the interband 
pairing, for cyx < c^2 Ucal , the minimum due to the nega- 
tive interference effect in 2Ai/T c at (fi — E e d ge )/^a = — 1, 
associated with the Feshbach-like shape resonance at the 
band edge, moves toward (/i — E e d ge )/<-^o = 0. The ratio 
2Ai/T c shows a maximum at (fi — E e d ge )/u>o = [dotted 
line in panel (c) of Fig. 8] for the critical value of the 
interband pairing, c c ^ tlcal ', 



IV. THE CASE OF THE SUPERLATTICE OF 
ATOMIC BORON LAYERS: DOPED MgB 2 

The present shape resonance superconductivity sce- 
nario in a multiband superconductor tuning the chem- 
ical potential near a band edge can be tested in the par- 
ticular experimental case of doped magnesium diboride. 
Doped MgB2 is a practical realization of a heterostruc- 
ture at atomic limit 16 being made of first units, super- 
conducting atomic boron layers (similar to graphcnc lay- 
ers), intercalated with second units (playing the role of 
spacers), hexagonal Mg layers, forming a superlattice of 
quantum wells with a period of 0.35 nm. 27,52 The elec- 
tronic structure near the Fermi level resulting from the 
material architecture is made of a first large Fermi surface 
(7r-band) and second small Fermi surface (cr-band) with 
different symmetry and spatial locations. The DOS at 
the Fermi energy in MgB2 is 0.12 states/(eV-atom-spin), 
where about one half (exactly 44%) of this value comes 
from the cr-bonds and the other half (exactly 56%) comes 
from the 7r-bonds. The chemical potential Ep in MgB2 is 
separated by about 700 meV from the band edge energy 
at the top Er of the cr-band and it is separated by about 
300 meV from the energy of the 3D-2D ETT at the E A 
point in the band structure. 

The chemical potential can be tuned near the band 
edge of the sigma band (at the critical point E e d ge = Er 
in the band structure) and across the 3D-2D ETT (at 
the critical point E3D-2D = E A hi the band structure) 
by chemical substitution of Al 27 ' 53 ~ 55 and Sc 56 for Mg or 
C for B. 57 ' 58 The different superconducting gaps in both 
bands and the critical temperature have been measured 
as a function of the chemical potential across the 3D-2D 
ETT by several groups in Al doped samples, 59 ~ 63 and in 
carbon doped samples. 6 The physical properties around 
the 3D-2D ETT are plotted here for the case of doped 
MgB2 as a function of the reduced Lifshitz parameter 
V = (E r - E F )/(E A - Ep), where £o=E A - E T is the 
energy dispersion in the c-axis direction due to electron 
hopping between the boron layers (400 meV in MgB2 and 
it changes with chemical substitution x). The influence of 
the proximity to an ETT 65 on the anomalous electronic 
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FIG. 9. (Color online). Panel (a). The supercon- 
ducting critical temperature in Al for Mg substitution 
(Mgi- x Al x )B2,(open circles) 53 , in Sc for Mg substitution 
(Mgi- x Sc x )B2,(open triangles) 56 and for the C for B substi- 
tution in the Mg(C x Bi- x )2 system (open squares) 58 as func- 
tion of the Lifshitz parameter "z" = {Er — Ef)/{Ea — Er) 
for diborides. 53 ' 56,58 Panel (b). The gaps Ai = A CT (open 
symbols) and A2 = A„- (filled symbols) in Al for Mg sub- 
stitution (Mgi- x Al x )B2, 59 and for the C for B substitu- 
tion in the Mg(C x Bi- x )2 system 64 as function of the Lif- 
shitz parameter "z" = (Er — Ef)/(Ea — Er) calculated for 
aluminum 53,54 and carbon 57 substitutions based on theoret- 
ical calculations and structural data. Panel (c). The BCS 
ratios curves 2Ai/T c (solid line) and —2A2/T C (dashed line), 
as functions of T c , extracted from the experiments measuring 
in the same run the two gaps well below the critical tempera- 
ture and T c . 59 ~ 63 Carbon doped samples, 64 (open circles). Al 
doped samples from Daghero et al. 59 (open squares), and from 
Samuely et al. 61 (triangles). The experimental data are com- 
pared with the calculated curves in the case C22/C11 = 2.17, 
C12/C11 = —0.43, where — C12 = 0.10 is lower than the value 
(— c c & tlcal ) = 0.52, where the critical temperature, and the 
gaps are calculated as a function of (p — E2,l)/£o- 



and lattice properties of MgB2 is shown by the anomalous 
pressure dependence 66 and Raman lineshape 67 of the E 2 g 
phonon mode. 

The multiband superconducting phase of MgB2 in the 
clean limit has been studied theoretically in detail 47, 68 ~ 70 
and its properties have been extracted from experimen- 
tal data. 64,71 The ratio between the intraband electron- 
phonon (p-h) coupling in the cr-band, C22, and in the 
7r-band, en, is in the range 2 < C22/C11 < 2.6 and the 
ratios between the interband coupling constant and the 
intraband coupling constants are 0.15 < C12/C22 < 0.3 
and 0.4 < C12/C11 < 1. The system is therefore in 
the regime of strong coupling in the second band and 
weak interband coupling discussed above. The inter- 
band coupling is smaller than the critical interband cou- 



pling 



^critical 
'12 



Y / ciiC22/0.656: from the published 



data we find 0.2 < c\ 2 j{- 



„critical 
-12 



) < 0.5. Therefore 



the shape resonance regime for the doped MgB2 system 
is well described by the typical case in Fig. 6 shown 
as the curve 5 where C22/C11 = 2.17, C12/C11 = —0.43 



and Ci 2 /(- 



„critical 



) = 0.19. For diborides there is a 



lack of a focused experimental test for the so called s± 
pairing, driven by the repulsive Coulomb interaction in 
the two band superconductivity with negative interband 
coupling. This will lead to condensate wavefunctions 
with opposite signs as it is proposed here. Therefore 
we hope that an experiment like for pnictides 49 will be 
run in the near future. However the critical tempera- 
ture and the absolute values of the gaps calculated here 
does not depend in the sign of the interband coupling. 
The critical temperature (panel a) superconducting gaps 
(panel b) in doped MgB2 are shown in Fig. 9 as a 
function of the Lifshitz parameter " z" . The Lifshitz pa- 
rameter "z" has been calculated at each chemical dop- 
ing value x in the case of the Al substitution for Mg 
{Mg 1 _ x Al x )B 2 , 53 ~ 55 for the case of Sc substitution for 
Mg in the (Mg 1 _ x Sc x )B 2 system 56 and for the C for B 
substitution in the Mg(C x Bi_ x ) 2 system 57,58 based on 
the measured variation of the lattice parameters and the 
charge density using band structure calculations. 53 ~ 55 > 57 
Fig. 9a shows the universal scaling of the critical tem- 
perature T c , as a function of the Lifshitz parameter 
"z" while the doping mechanism is quite different for 
the aluminum, 53 ' 54 scandium 56 and carbon 57 ' 58 chemi- 
cal substitutions in diborides. It is clear that the critical 
temperature curves collapse on the same curve for all 
cases when they are plotted as a function of " z" . The 
experimental data are plotted in the same figure with 
the theoretical calculations for the generic T c and gaps 
behavior for the case 5 in Fig. 6 and it is possible to 
appreciate the good agreement between theory and ex- 
perimental data. In Fig. 9c, the BCS gap ratios 2Ai/T c 
and — 2A 2 /T C of aluminum doped and carbon doped di- 
borides are plotted as functions of T c and are compared 
with the same theoretical case discussed above. The 
agreement between our theoretical calculations and the 
available experimental data is satisfactory for the depen- 
dence of T c , the gaps, and the BCS ratio in a wide range 



11 



of aluminum doping of magnesium diboridcs. The differ- 
ence between the experimental data and the present cal- 
culations is related with the fact the intraband coupling 
in the second band is likely to depend on the chemical 
potential, 67 whereas we assume it to be constant, to re- 
duce the number of parameters. Further work is under 
progress for a quantitative investigation of shape reso- 
nances in doped diborides. These results show clearly 
that superconductivity in doped diboride superlattices 
remain in the clean limit for interband pairing although 
the large number density of impurity centers in doped 
samples and the structural phase separation. 67 Indeed, it 
has been pointed out that the shape resonance mecha- 
nism can be favored by phase separation and especially 
for fractal distributions in supcrstripes. 21-23 



V. CONCLUSIONS 

The theoretical analysis discussed in this paper was 
motivated by the crucial observation that interband pair- 
ing may yield a sizable enhancement of the critical tem- 
perature in multiband materials, displaying superconduc- 
tivity in the clean limit, possibly leading to room tem- 
perature superconductivity, 16 and is the first step toward 
the understanding of the role of shape resonances in the 
superconducting gaps in multilayer systems. Although 
limited to the BCS mean-field approximation, our ap- 
proach captures the main ingredients of the physics of 
these systems. In particular, once the relevant details 
of the wave functions and of the band dispersions of the 
multilayer structure are introduced in the coupled BCS- 
like equations for the gaps in the various bands and for 
the conservation of the electron density, the resonance 
physics can be investigated in the very same spirit of the 
BCS-BEC crossover in fermion gases, where a similar ap- 
proach proved to be predictive. 

In real materials, the shape resonances can be con- 
trolled with material science techniques, tuning the pair- 
ing strength and tailoring the thickness of the spacer 



layers, to adjust the electron hopping between supercon- 
ducting layers. In our theoretical approach, the crucial 
parameter to take into account the effects of hopping vari- 
ations on the characteristics of the shape resonance, is the 
dispersion £. Resonance effects are emphasized when £ is 
on the order of the energy cutoff ujq of the effective pair- 
ing interaction leading to Cooper pairing. In this paper, 
we devoted our analysis to the case £ = loq , but of course 
£ can be varied around the optimal amplification value, 
and an analysis of the dependence of resonance effects on 
the value of £ will be the object of future investigation. 

The main outcomes of our theoretical scenario for su- 
perconducting multilayers or multigap superconductors 
have been tested here on aluminum-doped magnesium 
diboride superconductors. The quantitative comparison 
between our theoretical results and available experimen- 
tal data for this material is satisfactory and permits to 
identify, e.g., the relevant range of values for the pairing 
coupling constants. 

Our approach can be easily extended to pursue a pre- 
liminary investigation of the superconducting properties 
of other multiband system where the chemical potential 
is tuned near a ETT in one of the bands. For instance, it 
could be applied to Fe-based 1111, 122 and 11 pnictides. 

On the basis of our analysis, we may start to indicate 
a possible roadmap for the discovery of novel HTS. 16 
The conditions for higher temperature superconductiv- 
ity, or even room temperature superconductivity, might 
be met in graphene bilayers and graphane, 72 ' 73 where 
the exchange-like interband interaction could be even 
stronger than the intraband Cooper pairing, as also pre- 
dicted for superlattices of carbon nanotubes. 22 
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